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ABSTRACT 

Loop Heat Pipes (LHPs) have proven themselves as 
reliable and robust heat transport devices for spacecraft 
thermal control systems. So far, the LHPs in earth -orbit 
satellites perform very well as expected. Conventional 
LHPs usually consist of a single capillary pump for heat 
acquisition and a single condenser for heat rejection. 
Multiple pump/multiple condenser LHPs have shown to 
function very well in ground testing. Nevertheless, die 
test results of a dual pump/condenser LHP also revealed 
that the dual LHP behaved in a complicated manna- due 
to the interaction between the pumps and condensers. 
Thus it is redundant to say that more research is needed 
before they are ready for 0-g deployment. One research 
area that perhaps compels immediate attention is the 
analytical modeling of LHPs, particularly the transient 
phenomena. Modeling a single pump/single condenser 
LHP is difficult enough. Only a handful of computer 
codes are available for both steady state and transient 
simulations of conventional LHPs. No previous effort 
was made to develop an analytical model (or even a 
complete theory) to predict the operational behavior of 
the multiple pump/multiple condenser LHP systems. 
The current research project offered a basic theory of 
the multiple pump/multiple condenser LHP operation. 
From it, a computer code was developed to predict the 
LHP saturation temperature in accordance with the 
system operating and environmental conditions. 

INTRODUCTION 

Most of the analytical modeling efforts focused on the 
temperature prediction of single pump/single condenser 
(conventional) LHPs after a steady state condition was 
reached [1, 2, 3], In these steady state models, energy 
balance among the heat leak across the primary wick, 
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return liquid subcooling, and environmental parasitics 
was generally utilized to determine the loop saturation 
temperature for a given set of operating and ambient 
conditions. Predictions made by tbe models, for the 
most part, correlated well with the corresponding test 
data. 

A transient model for conventional LHPs was recently 
developed under the direction of NASA Goddard Space 
Flight Cotta- [4], The governing equations were, of 
course, derived from the basic conservation laws of 
mass, momentum and energy. However, observations 
of unique LHP operational characteristics during testing 
also assisted the formulation of the LHP theory. The 
transient LHP model (FORTRAN computer code) was 
independently verified by several users with different 
LHPs [5, 6, 7]. The predictions-versus-data comparison 
was outstanding (error bounds < 2°C) even under severe 
operating conditions (e.g. rapid change of power/sink). 
But perhaps more important than that, the computer 
runtime was extremely efficient taking less than 2 
seconds for a 1.8GHz Pentium IV laptop computer to 
simulate 10 hours of LHP operation. It did not exhibit 
any kind of numerical instability, either. 

In the follow-on research, it was therefore logical to 
extend the aforementioned transient LHP theory/model 
for the analytical simulation of the NASA/GSFC dual 
pump/condenser LHP system. As it turned out, the 
methodology in formulating the theory and deriving the 
governing equations was more or less straightforward 
(albeit tedious). But the mechanics of solving tire 
equations became so complicated that a new computer 
code had to be written from scratch for the multiple 
pump/condenser LHPs. Two versions of the code were 
developed in the current project. The steady state 
version was completed and verified with available test 
data. It is therefore the main focus of this papa. The 
theory, governing equations, solution scheme, and data 
correlation will be presented. The transient version was 
also finished and currently in the process of data 
comparison and evaluation. The results will be 
presented at a public forum in the near future. 
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NOMENCLATURE 

A® * cross-sectional area of condenser i. 
c p = liquid specific heat of working fluid. 

D® = inner diameter of condenser i. 

D® w k wick inner diameter of pump i. 

D® w = wick outer diameter of pump i. 

Gr-a s thermal conductance from reservoir i to 
ambient. 

h 2<p = film coefficient for two-phase pipe flow. 
h L m film coefficient for single-phase pipe flow, 
k^ s effective thermal conductivity of wick i. 
s length of condenser i. 

= two-phase length of condenser i. 

= length of wick i. 

m£f = mass flow rate in condenser i. 
m® = mass flow rate in pump i. 

N c = number of parallel condensers. 

N p = number of parallel capillary pumps. 

Q® = heat leak across wick i. 

Q®max 55 maximum heat leak across wick i. 

Q® = heat removal by condensation in condenser i. 
Q® s heat input to capillary pump i. 

Q c = heat removal by condensation at condenser. 
Qr-a = heat loss from reservoir i to ambient. 

Q®. s return liquid subcooling to pump i. 

Q sc, max = maximum return liq. subcooling to pump i. 
4m® 

Re == — S— = Reynolds number of fluid flow. 

"D®p 

T® = temperature of return liquid to pump i. 

T sat s saturation temperature of LHP. 

T sink s sink temperature. 

T„ = ambient temperature, 

x = vapor quality in two-phase flow. 

Greek 


v L = liquid kinematic viscosity of working fluid. 
p t = liquid density of working fluid. 
p v e vapor density of working fluid. 

X = latent heat of vaporization of working fluid. 
<D l = two-phase function. 

AP C == pressure drop across condensa's. 

AP® l = pressure drop in liquid portion of condenser i. 
AP® 2<1 = 2 -phase pressure drop of condenser i. 

AP®^ = pressure drop across flow regulator i. 

AP ® s pressure drop across meniscus of wide i. 
AT$ = temperature difference across wick i. 



slope of Tsat versus P S at curve. 


LHP GOVERNING EQUATIONS 

The derivation of the governing equations for a multiple 
pump/condenser is broken into two independent parts. 
One is for the heat rejection and fluid flow distribution 
among the parallel condensers. Note that capillary flow 
regulators, placed downstream of the condensers, may 
be needed to prevent vapor blow-through in case of a 
severe flow imbalance. The other set of equations i$ for 
the energy balance in the parallel capillary pump/hydro- 
accumulator assemblies. 

PARALLEL CONDENSERS 

The equations for steady state heat rejection/fluid flow 
distribution among the LHP parallel condensers were 
derived from three following basic principles: 

• pressure drop across each condenser is identical. 

• mass flow rate in a condenser is exactly equal to 
the condensation rate of heat removal. 

• when a condenser is folly open, the flow regulator 
provides additional pressure drop to make up the 
difference between foe total condenser pressure 
drop and its own frictional losses. 

As depicted in Figure 1, foe frictional pressure drop in 
each condenser is foe sum of pressure drops in foe two- 
phase and liquid-phase portions of the condenser, i.e. 


a s local vapor volume fraction of two-phase flow. (1)-.. AP C = AP ® 2q + AP ® L + AP^ fori=ltoNc 

s = porosity of primary wick. 

p L = liquid viscosity of working fluid. where both AP^ and AP® L are functions of foe mass 

p v = vapor viscosity of working fluid. flowrate m® and two-phase length L® 2(p of condenser 
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Figure 1 - Flow Distribution Among Parallel Condensers 


i. Darcy’s law for single-phase pipe flows is employed From Eqs. (1), (la), (lb) and (2), AP C is an explicit 

to calculate AP® L as follows: function of m® , i.e. 


(la)... AP® l 


L^-L®,^ 1 


fm®^ 


D® 


2p 


V A c ) 


for i = 1 to N c 


where, 

f 64 / Re laminarflow 

{ 0.316/ Re 1/4 turbulent flow 


For the two-phase pressure drop AP®^ , different two- 
phase models, depending upon the flow regime, can be 
used [8]. In each of these two-phase models, the heat 
transfer coefficient and pressure drop were correlated 
with respect to the local vapor void fraction, i.e. 




T > 

J f dx^ 

Idz J 


dx 


where. 


l _ Tsat ~ 
d z) m 


d z J m®X 


-JlD®h 


C “29 


= 


1-a 


The two-phase length L^ 2(? is also a fiinction of the 
mass flow rate m® , i.e. 


(3) ... AP C = F i (mg ) ) for i = 1 to N c 

There are N c equations in (3) for N c +1 unknowns: m® 
and AP C . To close the problem, employ the continuity 
equation for the mass flow at die condenser inlet: 

(4) ... ^m^mv 


Solving Eqs. (3) and (4) simultaneously will obtain the 
mass flow rates m® in the parallel condensers. Note 
dial when a condenser is folly open, vapor is prevented 
from flowing through that condenser by the capillary 
flow regulator. In other words, die flow regulator exerts 
an additional amount of capillary pressure to force the 
excess vapor to flow to other condensers as illustrated 
in Figure lb. 

PARALLEL CAPILLARY PUMPS 

Referring to Figure 2, for a multiple parallel pump LHP 
to reach steady state, the following conditions must be 
satisfied: 

( 5 )... 

JQ^max = Qsc.max + Qr-a fori = K 
IQ®max^Q^max+Qr-a fori = ltoN p andi*K 


( 2 )- 


m®A. 


L « = f— dx 

*D (i) fr -T Wh 

7LL, C V 1 SAT *SINK/0 “29 


Eq. (5) stipulates that the worst-case energy balance 
must be achieved for at least one of the pumps/hydro- 
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Figure 2 - Two-Phase Flow in Capillary Pump Core 

accumulators (i.e. pump K). The worst-case heat leak 
condition occurs when both the hydro-accumulator and 
pump core contain both liquid and vapor. The heat leak 
Q?max the amount of return subcooling QscImax 
in pump K are: 


( 6 )... 

Q 



AT® 


with C, = 


mjPc 


2nk®L® 


(7)- Q^,max =m? ) c P (T SAT -T®) 


Substituting Eqs. (6) and (7) into Eq. (5), will obtain: 



Figure 3 - Heat Leak Across LHP Wick 


For other pumps, the energy balance is still met under a 
different condition. Starting with Eq. (5), when the heat 
leak Q®max is less than the combined cooling effect of 
Q sc, max Q® A , then the vapor void fraction in the 

pump core and hydro-accumulator decreases, reducing 
both the heat leak Q® and die return subcooling Q®. 
as depicted in Figure 3. Note that the reduction of Q^. 
is due to the incomplete heat exchange between the 
return subcooled liquid in the bayonet tube and the 
pump core vapor (i.e. the liquid exiting the bayonet is 
still subcooled). The vapor void fraction continues to 
decrease in pump i until it readies a level at which the 
energy balance prevails. In other words, a two-phase 
condition may exist in the core and hydro-accumulator 
of pump i. However there is also a strong possibility 
that the pump core/hydro-accum ulator is filled with 
liquid if the return subcooling overwhelms the heat leak 
in pump i. 

In die steady state model, die worst-case heat leak is 
calculated to determine the LHP operating temperature 
T« t that satisfies the energy balance of each pump, i.e. 


(8)... Q"max-Q£,max 


+ Q£a 


pCO 


To achieve the Systran thermal equilibrium (or to satisfy 
Eq. (5)), the LHP operating temperature is die highest 
of all T<j> T . Or, 


m®c P 


(5a)... (D» w /D®)«-1 
where. 


AT®=m®c,(T SAT -T®)+ 

G&(Tsa T -TJ 


AT®> 


i 


5T' 
SP, 


AP® 


SAT 


The total pressure drop across the pump wick AP® 
depends strongly on the location of liquid-vapor front in 
the condenser, which in ton varies with the saturation 
temperature. Moreover, the thermo-physical properties 
of the working fluid are functions of temperature. Thus 
Eq. (5a) is usually solved iteratively to obtain the LHP 
operating temperature T S at- 


(9)~. T sat =max(l£> T ) for i=ltoN p 

STEADY STATE COMPUTER COPE 


A computer code, written in Visual BASIC, is used in 
conjunction with Microsoft Excel to predict die steady 
state operating temperature T sat of a multiple capillary 
pump/condenser LHP. For a given operating condition, 
die model goes through an iteration process to calculate 
T sat that satisfies Eq. (5). The iteration steps are 
described below: 

(i) guess loop saturation temperature T S 'JF SS . 

(ii) calculate thermo-physical properties (e.g. density, 

viscosity) of working fluid at . 
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but shared common transport lines. Each evaporator 
pump had its own concentric hydro-accumulator. The 
evaporator casings were made of aluminum tubing 
having 15.8mm (5/8”) O.D. and 76.2mm (3”) length. 
The hydro-accumulators (14.8mm O.D. x 81.8mm L), 
machined from 316L stainless steel tubing, were 
welded to the evaporators by way of bi-metallic joints. 
One evaporator had a titanium wick and the other had a 
nickel wick. Properties of the wicks are summarized in 
Table 1. 


The loop saturation temperature is therefore regulated 
(or controlled) by the hydro-accumulator that satisfies 
the Equation (5) identically. As far as the steady state 
modeling is concerned, the exact thermodynamic states 
of the other hydro-accumulators are not required and 
therefore omitted in die current calculation. Note that 
the same cannot be said about the transient version of 
die LHP model in which die thermodynamic conditions 
of all hydro-accumulator had to be computed at every 
time step. 

MODEL VERIFICATION 

The verification of the LHP steady state computer code 
was made by correlating die model predictions versus 
available test data from a dual pump/condenser LHP as 
shown in Figure 4. The dual LHP test article, built by 
the Dynatherm Corporation, consisted of two parallel 
capillary evaporator pumps and two parallel condensers 


Table 1 - Wick Properties of NASA/GSFC Dual LHP 


Property 

Pore Radius 
Porosity 
Permeability 
Thermal Conductivity 


Titanium Wick Nickel Wick 
3.0microns 0.6microns 

45% 60% 

7.4xl0' 14 m 2 1.4xl0' 14 m 2 

7.2W/m-K 5.8W/m-K 


The vapor and liquid lines were made of stainless steel 
tubing, having the respective outer diameters of 2.2mm 
(3/32”) and 1.6mm (1/16”). The length of both vapor 
and liquid lines measured approximately 1,168mm 
(46”). Each condenser was formed by sandwiching a 
serpentine 2.2mm O.D. x 1,016mm L stainless steel line 
between two 82.6mm x 177.8mm (3.25”x7”) aluminum 
plates. A flow regulator made of 20pm-pore capillary 
screens was placed downstream of the condensers. The 
flow regulator prevented vapor from returning to the 
evaporator pumps when one of the condensers was fully 


(iii) solve Eq. (3) to determine heat removal/fluid flow 
distribution among parallel condensers. 

(iv) calculate pressure drops in all LHP components to 
obtain total pressure drops across each pump wick. 

(v) compute Q®^, Q% MAX from Eqs. (6) and (7) 
and Q® A , then determine if the energy balance in 
pump i satisfies Equation/Inequality (5). If it does, 
exit iteration loop. Otherwise, go back to step (i) 
to Tepeat process until it does. 






NASA/GSFC Dual I-HP Operating Temperature 




280 1 1 1 1 1 1 1 1 1 1 

0 10 20 30 *0 50 60 70 80 90 

Power Input Per Evaporator (W) 


Figure 5 - Model Predictions vs. Test Data 

utilized, and hence served to balance the fluid flow 
between the two condensers. The loop was charged 
with 15.5 grams of anhydrous ammonia. 

Electrical heaters were attached to both evaporators and 
hydro-accumulators and the heater power inputs were 
regulated by separate power supplies. The two 
condensers were clamped to two cold plates; each 
cooled by a chiller. Sixty thermocouples were used to 
monitor the loop temperatures. Notice that many 
thermocouples were installed on the liquid line between 
the two hydro-accumulators in order to detect any 
interaction between the two during last transients. A 
data acquisition system consisting of a data logger, a 
personal computer, a CRT monitor, and Lab view 
software was used to monitor and store data. The data 
was updated on the monitor and stored in the computer 
every second. 

Figure 5 presents the temperature comparison between 
the test data and model predictions for the case in which 
the same heater power was applied to each evaporator 
pump. In the tests, two sink conditions (-10°C and 
+10°C) were imposed on the loop operation. The 
correlation was excellent. The opmating temperature 
from the model simulation came within ±2°C of the test 
data in this test series when the power input was less 
than 60W/pump for the +10°C sink and 50W/pump for 
the -10°C sink. For power inputs higher than these 
levels, the model showed that the system pressure drop 
exceeded the capillary limit of the titanium wick. Note 
that this condition only implied that vapor vented into 
the pump core through some area of the wick (partial 
dryout) but might not necessarily lead to a complete 
dryout of the titanium wick (i.e. pump #1 deprime). 
Nevertheless the actual LHP saturation temperature was 
progressively higher than those predicted by the model 
as the power input increased beyond 60W/pump. In 
other words, even though pump #1 did not foil, the 


NASA/GSFC Dual LHP Operating Temperature 
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Figure 6 - Model Predictions vs. Test Data 

vapor leakage through the wick had the same effects as 
the heat leak Q® . Both test data and model predictions 
indicated that the hydro-accumulator #1 of the titanium 
wick pump actually controlled the loop saturation 
temperature. This made perfect sense, since the 
titanium wick has a higher effective thermal 
conductivity than that of the nickel wick, resulting in a 
higher amount of heat leak Q® into the hydro- 
accumulator #1. Under the same heat input, each pump 
drew about the same amount of return liquid subcooling 
Q®. . Thus according to Eq. (5), the hydro-accumulator 
#1 had to drive the LHP to a higher temperature to 
satisfy its energy balance. 

Figure 6 shows the test results and model predictions of 
die uneven power tests. During this series of tests, both 
sink (condenser) temperatures were set at 0°C. The 
power inputs to the evaporators were varied with a 25W 
increment but the combined power was maintained at 
100W. As expected, the evaporator that had less heat 
input controlled the loop saturation temperature (since 
it drew less subcooling Q§. back to the pump). The 
model prediction of the loop saturation temperature was 
also very good for the uneven power tests. In addition, 
it correctly predicted the hydro-accumulator that 
controlled the LHP temperature for a given condition. 

CONCLUSIONS 

Loop Heat Pipes have become increasingly popular 
with spacecraft thermal engineers for their operational 
reliability and robustness. Well over 100 conventional 
LHPs have been utilized for spacecraft Thermal Control 
Systems (TCS). Thus multiple evaporator/condenser 
LHPs are the next logical step in the LHP utilization. 
Indeed several flight projects are considering Ihem as 
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the primary TCS. Appropriately, the thermal engineers 
need accurate LHP models to simulate die thermal 
system performance particularly on-orbit behaviors. 

The LHP theory of conventional LHP operation was 
extended to multiple pump/condenser systems, and 
from which, a steady state computer code was written. 
The computer code was simple and easy to use. As 
presented in this paper, the computer predictions agreed 
well with available test data. Nevertheless, the 
modeling effort should continue to provide a better 
simulation tool for space applications especially in 
transient predictions. 
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